Hyperfine Effects 
and Large I/tuq Corrections 



ON 
On 



OO 
(N 

(N 
> 

(N 



B. Holdom* and M. Sutherland^ 

Department of Physics, University of Toronto 
60 St. George St., Toronto, Ontario, Canada M5S 1A7 



UTPT-92-24; [hep-ph/92 12277 



Abstract 



Effects associated with hyperfine mass splitting between pseudoscalar and 
fN| \ vector mesons may induce surprisingly large I/tuq corrections in the heavy 

quark expansion. We demonstrate this in a relativistic quark model by calcu- 
li . lating to all orders in l/mg, with and without hyperfine effects, and comparing 
to the first order results. Total corrections of 30% or more are quite possible 
O in the decay rates for B — > Dlv and B — > D*lv near zero recoil. 

6-;, 

The subject of semileptonic weak decays of mesons containing a heavy quark has 
recently attracted much attention. In the limit in which the heavy quark mass vhq 
becomes infinite the theoretical treatment of these decays is greatly simplified, 



potentially reducing the uncertainty in the extraction of the Kobayashi-Maskawa 
elements. But more insight into the magnitude of the corrections to this limit is 
required. 

We present here a relativistic quark model which provides a representation of 
the decay amplitudes to all orders in 1/mq. In we show that the model correctly 
incorporates all heavy quark symmetry relations at zeroth and first order in the 1 / vtlq 
expansion. We will use the model to study a source of corrections to the heavy quark 
limit which has not yet been seriously considered. This is the spin symmetry breaking 
effects due to gluon exchange between the heavy and light quarks. These hyperfine 
effects split the vector and pseudoscalar masses, but they will also cause a distortion 
of the light quark wave function in a way which depends on the spin of the heavy 
quark. 
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Our relativistic model does not deal with gluon exchange directly, but it relates 
the distortions in light quark "wave functions" to the meson mass splitting. In || 
we explored the sensitivity of amplitudes describing B — > Dlv and B — > D*lu to this 
hyperfine splitting at first order in the heavy quark expansion. Hyperfine splitting 
had no effect on any of the amplitudes at zero recoil. This gives some motivation for 
the consideration of 0(1/ mq) effects, in addition to the fact that the main source of 
corrections affecting the extraction of V c b occurs at this order. 

A systematic analysis of 0(l/m 2 Q) corrections in the heavy quark effective theory 
has recently been made in [|J. Some attempt was also made to estimate the quantities 
relevant to semileptonic decay at zero recoil, and to do this the authors appealed to 
the ISGW nonrelativistic quark model. 0] But the ISGW meson wave functions are 
"spin averaged" ; hyperfine effects are not included. The conclusion in || is therefore 
consistent with ours: if hyperfine effects are ignored then 0(1/ tuq) corrections to 
semileptonic decay rates at zero recoil are small. But we are led to a very different 
result when hyperfine effects are included. 

Our model Lagrangian will contain terms which couple the heavy meson fields to 
the heavy and light quark fields. 

-Q(x){i lh P(x)F P (x -y)+ 1 ^(x)F v (x - y)}q(y) + h.c. (1) 

Q and q are the heavy and light quarks contained in the heavy pseudoscalar and 
vector mesons P and V. Fpy(x — y) are damping factors which will suppress large 
momentum flow into the light quark. This is the essential physical effect of the light 
quark wave function. In momentum space we take 

where k is the momentum of the light quark. The quantities Z P y and A P y are com- 
pletely determined for each meson in a manner to be described below. In particular 
the fact that A P is not equal to A v reflects the spin dependent distortion of "wave 
functions" mentioned above, and Apy will be fixed by the physical, spin dependent 
meson masses. Such a connection between the distortion and the meson mass split- 
ting will be present in any quark model of heavy mesons, whether or not the model 
explicitly relates hyperfine effects to gluon exchange. 

In we considered the effect of raising Fpy in (|2]) to powers n other than unity. 
We found that meson decay constants were the physical quantities most sensitive 
to the high energy behavior of the vertex factors, and we found that they decrease 
for increasing n. To obtain realistic meson decay constants, values of n close to but 
greater than 1 are preferred. Although n — 1 yields a divergent vector meson decay 
constant, the n — 1 case modified by inverse powers of logarithms would probably 
also be acceptable. Such modifications would have minimal impact on the quantities 
of interest in this paper, and we shall focus on the pure n = 1 case as given by (0). 
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The rest of our model Lagrangian consists of standard kinetic and mass terms for 
the quarks. The heavy mesons are composed of quarks, and thus to avoid double 
counting of dynamical degrees of freedom the model Lagrangian does not contain 
meson kinetic terms. Neither does it contain mass terms or meson self-couplings; 
all these terms are to be modeled via the quark loop. Various other fields may 
be coupled to the quarks, depending on the application. For example to describe 
semileptonic decay we will couple the heavy quarks to the appropriate external gauge 
fields. Processes involving pions or kaons may be described by coupling these fields 
to the light quarks. 

The quark masses are the only parameters of the model. For the heavy quarks 
we use the values m& = 4.8 GeV and m c = 1.44 GeV. For the light quark we use 
its constituent mass m q ; but here there is some uncertainty. Since we expect that 
the actual momentum dependent constituent quark mass has fallen somewhat at 
momenta typical in the loop, we take a (momentum independent) m q « 250 MeV. 
We will also explore the sensitivity of results to changes in m q . 

The free quarks may be integrated out of the theory and this yields an effective 
theory of heavy mesons. The form of such a theory has been well described in the 
literature. |5| This is a theory of heavy mesons propagating on mass shell or close to 
mass shell. It does not describe diagrams with closed loops of heavy mesons. It may 
describe pion loop corrections with an internal heavy meson line since the effective 
cutoff on such a loop is far below the heavy meson mass. 

Here we note a peculiar feature of our model. It turns out that upon carrying out 
the procedure just described, the resulting effective heavy meson Lagrangian has an 
overall nonstandard minus sign. But this sign, which would be disastrous for a normal 
quantum field theory, turns out to be an unphysical sign for this effective theory. This 
is due to the absence of heavy meson loops. With each heavy meson vertex and each 
heavy meson propagator having an additional minus sign, the result is that a diagram 
receives a sign (— 1) J where j is the number of heavy meson lines running through 
the diagram. Such a sign is of no consequence for physical quantities calculated from 
these amplitudes.^] The amplitudes of interest here are determined by the tree level 
two-meson-external-gauge-field vertices with heavy mesons on shell. In the following 
we will simply omit the minus sign to coincide with standard conventions. 

We denote the proper pseudoscalar and vector self-energy graphs by iY P {p 2 ) and 
— ig^Tvip 2 ) + • • •, where the ellipsis denotes p^p u terms. Because the model treats 
the quarks as free, the Tpy become complex above the threshold at p 2 = (mQ + m q ) 2 . 
We will ignore this consequence of free quarks and drop all such imaginary pieces of 
quark loop graphs.0 The meson masses are defined by Y py(m 2 PV ) = 0. The mass 

^^One may wonder about the existence of negative energies. But terms of higher order in deriva- 
tives and heavy quark mass would be relevant for the Hamiltonian. 

2 The imaginary parts may in fact be kept for all quantities we calculate; the result would be a 
consistent description of mesons coupling to free quarks. Our dropping the imaginary parts in the 
end is an effort to incorporate confinement. In this connection we note that the confinement physics 
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nip Ap Zp 


V 


m v Ay Z v 


B 
D 


5.279 0.6214 1.1026 
1.869 0.5367 0.7446 


B* 
D* 


5.325 0.7015 1.2896 
2.010 0.7623 1.1927 


B 
D 


5.31 0.6661 1.1920 
1.97 0.6739 0.9808 


B* 
D* 


5.31 0.6789 1.2432 
1.97 0.7045 1.0878 



Table 1: Input meson masses and resulting values of Apy and Zpy with (top half of 
table) and without (bottom half) hyperfine splitting (in GeV). 

functions V B B * are plotted versus y/p 1 in Fig. 1. The zeros lie between the threshold 
kink and singularities at p 2 = (itlq + Apy) 2 . (We therefore obtain sensible meson 
masses only if Apy > m q .) The Apy are determined by requiring that the zeros of 
the meson mass functions coincide with the experimental meson masses. 

The constants Zpy are determined by normalizing in either of two equivalent 
ways. The Ward identities relate the slopes of the respective mass functions at their 
zeros to the q 2 = matrix elements of the heavy quark vector current between on-shell 
mesons. Normalizing these matrix elements to unity is equivalent to the condition 
r py (mpy) = 1, where the prime denotes differentiation with respect to p 2 . Thus, 

T Py {p 2 )^p 2 -m 2 P y (3) 

in the neighborhood of p 2 = m 2 PV , corresponding to standard meson kinetic terms. 
For the B and D sectors, the values of the input meson masses and the resulting 
values of Apy and Zpy are shown in the top half of Table |I| 

We are interested in the effects associated with hyperfine mass splitting, my—rap. 
In the model this originates largely in the different values of Ay and Ap. (There is 
also a little mass splitting even in the case Ay = Ap, with the vector mass slightly 
less than the pseudoscalar mass.) We will be comparing our results based on a fit 
to physical meson masses with results in the case of no hyperfine splitting. We will 
take the latter case to correspond to the values of Ap and Ay which yield both the 
pseudoscalar and vector masses equal to \mp + \my. The amounts by which the 
physical masses are shifted from this common mass is characteristic of the effect 
of hyperfine mass splitting. We display for comparison purposes the various input 
parameters for the case of no hyperfine splitting in the bottom half of Table [3]. 

We are also interested in comparing to results obtained at first order in the heavy 
quark expansion. || The heavy quark limit of the model is described by the conditions 

Apy -> A , A/m Q , A/m q fixed. (4) 

of QCD occurs at a lower momentum scale than the scale characterizing our model, Apy ~ 650 — 700 
MeV. 
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This common A is the same for all heavy quark flavors. At zeroth order in an expan- 
sion in A/rriQ, quantities depend only on the ratio A/m q and are independent of the 
way in which Apy — ► A. At first order in the expansion we write 



A P y = A[l-(5 P yh + g) (5) 




where 5p = 3 and Sy = — 1. In @ we study the dependence of first order corrections 
on the parameters A/m q , g and /i.Q 

The parameter h is responsible for most of the hyperfine mass splitting. We may 
choose a set of parameters in the first order model which produces an optimal meson 
mass spectrum. For m q = 250 MeV such a set is A = 667 MeV, g = —0.13 and 
h = 0.19 and it yields (B,B ) and (D,D*) masses within 0.2% and 2% respectively 
of the physical masses. 

The form factors for the semileptonic decays B — > D and B — > D* are defined as 
follows. 

(P 2 {v 2 ) K\ Pi(^i)) = ^M Pl Mp 2 [h + {uo){ Vl + v 2 )^ + h-(u)(v x - v 2 ) f ] (6) 
(V 2 (v 2 ) \V,\ P 1 (v 1 )) = ^MftM^M^W^X (7) 

(V 2 (v 2 )\A l ,\P 1 (v 1 ))= (8) 
-i^M Pl M v , 2 (u + l)h Al (uj)e* 2lx - (h A2 (u)v lfl + h M {uj)v 2 ^) £* 2 ■ v 1 

The v's are meson velocities, and uo = V\ ■ v 2 . 
At zeroth order in the heavy quark expansion, 

h+(u) = h v (uj) = h Al (u) = h A3 {uj) = £(u) , h-(u) = h A2 (u) = 0. (9) 

is the Isgur-Wise function. The first order corrections to the hi(u) require the 
introduction of only four additional universal functions. || In || we show that this is 
also true for our model, for any choice of parameters. This is a nontrivial test for any 
model of heavy mesons. 

We may calculate the hi(u) in the full model from the appropriate three point 
quark loop graphs. Our results are displayed in Fig. 2. For comparison we also 
display emerging at zeroth order in the model (using A = 667 MeV). We find 
that hv, A . A , Al ,+ {yj) display substantial shifts of varying amounts from although 
the shape of these functions is fairly similar. In addition /i_ iA2 {uj) deviate significantly 
from zero. 

We collect the full zero recoil values hi(l) in the first row of Table |2], and in the 
second row we collect the zero recoil values with no hyperfine splitting, h^ ohyp (l). 

3 Our definition of g and h here differs from that in B. 
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hi(l) 


1.401 


1.307 


1.155 


1.107 


-.138 


-.249 


h™ hyp {l) 


1.232 


1.149 


1.001 


0.994 


-.108 


-.274 




1.228 


1.152 


1 


1 


-.124 


-.270 


K(i) 


-2.21 


-1.95 


-1.71 


-1.83 


0.20 


0.51 


(/ir hyp )'(i) 


-1.76 


-1.54 


-1.32 


-1.53 


0.14 


0.51 


(^)'(i) 


-1.69 


-1.49 


-1.29 


-1.54 


0.13 


0.46 



Table 2: Zero recoil values of form factors hi (top half of table) and derivatives h\ 
(bottom half) with hyperfine splitting, without hyperfine splitting, and at first order. 



Most of the hi(l) deviate significantly further from their values in the heavy quark 
limit than do the /i, nohyp (l). Especially interesting are the quantities h + (l) and h^il) 
which are both still equal to unity in the heavy quark expansion at first order (Luke's 
Theorem |J). We see from Table |]that at higher orders, in the absence of hyperfine 
splitting, h+(l) and ^Ai(l) remain very close to unity. But the higher order correc- 
tions with hyperfine splitting cause significant deviation in these quantities away from 
the heavy quark limit. 

The zero recoil values from the first order model, hj st (l), are given in the third 
row of Table |[ In we found that these values are independent of g and h and 
thus independent of hyperfine splitting. The substantial difference between hi(l) and 
^nohyp^-Q . g therefore an 0(l/mg) effect. It is also interesting that the h} st (l) are 

very close to the /i° ohyp (l). This shows that were it not for the hyperfine effects, the 
corrections beyond first order would be very small. It appears that hyperfine effects 
completely dominate the 0(1/toq) corrections. 

In the bottom half of Table we give the values of the slopes of the form factors 
at zero recoil. By comparing the first two rows we see that the effect of hyperfine 
splitting is to make the slopes more negative. And the effect of hyperfine splitting at 
first order (that is, the shift in the (hj st )'(l) caused by nonzero h) is no more than 
about 1%. Thus again the main effect of hyperfine splitting occurs at 0(l/rrvn). 

In Table |]we show how the hi(l) change when the light quark mass m q is varied. 
Increasing it to 300 MeV increases the corrections even further. And the corrections 
remain substantial even if m q is lowered arbitrarily. 

We consider further the two quantities which are protected from first order cor- 
rections at zero recoil. The following relations are derived in [|3| . 
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m g (MeV) 


V 


^3 


A 1 


+ 




A 2 


300 


1.481 


1.390 


1.225 


1.156 


-.150 


-.259 


200 


1.363 


1.263 


1.123 


1.086 


-.132 


-.242 


10 


1.334 


1.193 


1.103 


1.082 


-.124 


-.208 



Table 3: Dependence of hi(l) on light quark mass m q . 

In a nonrelativistic quark model ^i(l) and ^2(1) are negative since they may be 
related to the deficit in a wave function overlap integral.^ This is not the case for 
our relativistic quark model. Our meson-g-Q vertex factors are analogous to wave 
functions, but they are inserted into a relativistic quark loop calculation. Our results 
explicitly show that as the two "wave functions" are distorted by hyperfine effects 
and made different from each other, the result is a positive correction. We have 
numerically isolated the 1/itlq corrections in (p70|,p7T|) and find A « (1.6 GeV) 2 , ^i(l) ~ 
(1.0 GeV) 2 , and #2(1) ~ (0.6 GeV) 2 . The total magnitude of the I/ttiq corrections in 
([l(JO) amounts to about 2/3 of the corrections to all orders. 

The A term in (|II]) may be written as 1.5A 2 /mfcm c ; the dimensionless coefficient 
here does not appear to be unreasonably large. In A is written as 

A= -A 1 + 2A 2 + ... (12) 

where the ellipsis represents corrections resulting from a double insertion of the 
chromo-magnetic operator. The A's are defined by 

Ai = — -(5m 2 p + 35m v ) and A 2 = — -(5m 2 p — 5m v ) (13) 



where 



m M - m Q = A + — + • ■ ■ 14 

ZlTLn 



We find in our model that A = 0.50 GeV, Ai = -(0.27 GeV) 2 and A 2 = (0.29 GeV) 2 . 
These A's are small contributions which largely cancel in A, and thus the bulk of 
contribution to A must originate in the double insertions of the chromo-magnetic 
operator. In || it was assumed that these contributions could be neglected. 

We now turn to the two processes which, by observing at zero recoil, have been 
proposed || as fairly model independent methods for extracting V^. In the differential 
decay rate at zero recoil for the process B — ► Dlv the following factor appears 



_ m B -m Dh _ 
m B + m D 



(15) 



4 We find from M £i(l) = l 2 {\) w -\rr? , i.e. 1/4 times those in 



7 



The deviation of this quantity from unity gives the correction to the heavy quark 
limit. This might be expected to be small since the only first order 1/rriQ correction 
occurs in h-(l) and this is suppressed by the prefactor. But inserting values from 
Table |2| gives a positive 38% correction with hyperfine splitting and a 9% correction 
without hyperfine splitting. For the decay B — > D lv the factor of interest is \\ia x (1) | 2 
which has no first order I/tuq corrections. Here the correction is 33% with hyperfine 
splitting, to be compared with essentially no correction for no hyperfine splitting. 

These numerical values for the corrections should be considered as illustrative. 
They depend not only on the light constituent mass, but also on the c and b quark 
masses we have chosen. And the corrections obviously reflect the simple form we 
have chosen for the meson vertex factors in (0). But the existence of significant 
corrections from hyperfine effects, with signs as given, seems to be generic to this 
class of relativistic quark models. 

In this paper we have argued by way of an explicit model that corrections to 
the heavy quark limit may be larger than previously thought. This is seen most 
graphically in Fig. 2 which depicts the form factors for semileptonic decay of heavy 
mesons. At zero recoil the larger than expected corrections appear at 0(1/ rug) in 
the heavy quark expansion. We have traced the origin of these corrections to the 
physics associated with hyperfine mass splitting. We hope that these results will 
encourage further study of hyperfine effects in nonrelativistic quark models and sum 
rule calculations. 
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Figure Captions 

Figure 1. Normalized B and B* mass functions for m q = 250 MeV. 
Figure 2. B — > £) and B ^ D* form factors /ii(a>) and Isgur-Wise function 
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Figure 2 



